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XXI-3. 
1. FORME NORMALI IN DINAMICA HAMILTONIANA 


Problema centrale nei sistemi dinamici hamiltoniani è l'analisi del 
la stabilità e della struttura geometrica delle orbite e può essere affrontato 
attraverso le forme normali e la teoria della coniugazione [1], [2]. 

Data una funzione H(x) definita su un compatto per, detto spa- 


zio delle fasi, le soluzioni delle equazioni di Hamilton 


O, I 
dx Ela 
(1.1) frode ica dI) ( ) 
dt dx 21 0 
dove I è la matrice identità dxd, definiscono una applicazione 35 RXD + D 
detta flusso hamiltoniano. I flussi hamiltoniani x' =2€.(x) hanno vincoli geo- 
metrici espressi dalla condizione detto di canonicità o simpletticità 


si gres: 
(1.2) MIM=I , Mx x 


dx 


Consideriamo un altro flusso x = F(X) definito da 


(1.3) = “ F(x,e) 


e notiamo che se x = (q,p), X = (Q,P) la forma differenziale 
(1.4) df = pdq + QdP 


è esatta e quindi il flusso è implicitamente definito da 


of of 
LE q=®= 


(1.5) p = 3q 5P 


Consideriamo quindi la trasformazione indotta da F. suH 


(1.5) H(#_(X)) = K(X,e) 
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si dice allora che H(x) e K(x)sono coniugati da 7. Se invece 
€ 


(1.6) H(F_(X)) = H(X) 

si dice che H è invariante. Se il flusso F, è autonomo cioè di = 0 allora 
e 

si ha anche 

(1.7) F(97,(x)) = F(x) 


ed i due flussi commutano 


bI 


- 


d 
Forme normali. Sia ora D = T x S dove S è una sfera in Ri e La 
il toro di dimensione d e sia x = ($,j) con set, jeS. Consideriamo inoltre 
jl gruppo G prodotto tensore delle rotazioni nei piani di fase. Se RE G con 


aeR 
(1.8) RX = (9 +0 .Î) 


Considereremo anche un sottogruppo 6, di G i cui elementi sono A con la re- 


strizione 


k= 0 belca dimL=r.. 
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Diciamo che H è in forma normale rispetto a G se 
È vRe 
(1.9) H(R x) H(x) Re G 


e si parla di forma normale non risonante. 


Diciamo che H è in forma normale rispetto a GL se vale 
(1.10) H( RA) = H(x) VA 6, 


e si parla di forma normale risonante. 


Se H è în forma normale non risonante allora 
(1.11) H= H(jjpccd,) 
ed il suo flusso A, è dato da 
(1.12) Hyx = (6 + wt, j) w= 77 
Se la frequenza w è non risonante cioè 
(1.13) uk#0  vkezio 


allora ogni orbita è densa su 1°; se invece w è risonante rispetto a Let 


con dim. L= r 


(1.14) w-ek=0 keL , w-k#0 VkéL 
< SIANO d-L 

allora ogni orbita è densa su un toro T°*. 

Osserviamo che se H è in forma normale rispetto a 6, può essere 


scritto nella forma 
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(1.15) og e E Hx(3) g1%*4 


dont 


keLl 


Consideriamo ora un hamiltoniano H che sia "vicino" ad un hamil- 


toniano Ho în forma normale non risonante 


(1.16) H(x,e) = H (5) +e Wo, i) 
e indichiamo con w = A la frequenza di Ho 


Le condizioni perché esista una trasformazione simplettica 
x = F(X) con x = (9,j) , X = (0,4) , è eîl che coniughi H con una forma normale 
H 


(1.17) H(#(X)) = H(X;e) 
sono espresse dal seguente 


Teorema. Se w è non risonante allora H è coniugabile con una forma 
normale (non risonante) rispetto al gruppo G; se w è risonante rispetto a L al 
Jora H è coniugabile con una forma normale (risonante) rispetto a GL: 

La condizione di coniugazione è verificata dalle serie formali in 


e se V è un polinomio trigometrico. 


Dimostrazione. Consideriamo una funzione generatrice F(9,J;e) per 
definire la trasformazione canonica F. e scriviamo il suo sviluppo in e nella 


forma 
(1.18) = F(6,0se) = 90 +D e" F,(6,0) 
nzl 


L'equazione di coniugazione diventa allora l'equazione di Hamilton-Jacobi e 


si scrive 
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,3F <A AIR, PRIN 
(1.19) Ho) +e V(ò, 30) H($,4;e) 
Sviluppando la 1.19 in serie nella variabile e e confrontando i 


ERIN . î . > 
termini di ordine e si ottiene îa ricorrenza 


aF 
(1.20) wo (0,9) + Q,(6,0) = (0,0) 


dove 


(1.21) = è (5,6) 


e Qh dipende solo da termini noti all'ordine n. 


(3), fi 


Se indichiamo con Fnk09)> Q nik 


(J) le componenti di Fou- 
nyk 


rier di F_(6,0), Q,(9,3) e A_(9,9) si trova 


(1.22) iw » k Frck * Mak = LR 


Se w è non risonante si ha 


H,=0  k=0 =Q F.=i Mk k#0 
n,k i 0) n,0 ? n,k we k 
mentre se w è risonante rispetto ad L si ha 
Bk =0 kéL, Uh “ak kel 
Q 
E n,k 
Foik i peri » KkÉEL 


Le componenti A K con k L rimangono arbitrarie. 
È] 


Mappe hamiltoniane. Chiamiamo mappa una applicazione simplettica 
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T ed orbita di T la successione n = Tx} dove 1° = ToTo..oT. Le mappe sono 


associate ai flussi mediante la sezione delle orbite del flusso con una super- 


; / [) 
A a+. 
7) nz CA 


Anche per le mappe si definiscono le forme normali rispetto al 


ficie (Poincaré) 


gruppo G del prodotto diretto dalle rotazioni nei piani coordinati 
(1.23) T(R (x)) =R(T(x)) 


Definendo XK = Vi, sin de PL = Vas, cos dk una mappa x' = T(x) in forma nor 


male può essere scritta nella forma 
(1.25) dr dti) i Fix k=1,...;N 


oppure posto Z = XK + iP, nella forma 


ia,(z3z}:-- 32 z*) 


ina nn - 
(1.26) ze zx k=1,...N 
pere 
dove 9,(4) = sj (j) 


Data dunque una mappa lineare perturbata 


= el “kz +f 


(1.27) zi K p(2:2*) 


ove la frequenza della parte lineare è costante, la forma normale con cui può 


essere coniugata è una rotazione i cui angoli dipendono dalle distanze dall'o 
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rigine in ciascun piano coordinato (ciò se w è non risonante) [3], [4]. Tali 
angoli di rotazione 2 sono dunque una caratteristica intrinseca di ciascuna 
orbita se n=2(j) è 1>1 ed è questa la causa delle difficoltà del problema 
della coniugazione in meccanica hamiltoniana e che rende la sua soluzione pos- 
sibile solo localmente cioè per ciascuna orbita. 

Un problema che le caratteristiche analoghe senza avere questo 
aspetto di località si incontra nella teoria della coniugazione per mappe olo 


morfe. 
2. CONIUGAZIONE PER MAPPE OLOMORFE 

CI limiteremo allo studio delle forme normali per mappe z' = F(z) 
dove F(z) è olomorfa in un disco |z| < R di C, ed ha nell'origine un punto fis 
so cioè F(0) = 0. 

SIa allora R(z) = Xz una trasformazione lineare. Diciamo che F_ è 
in forma normale rispetto al gruppo delle trasformazioni lineari che preser- 


vano l'origine se 


(2.1) FoR=RoF 


Se scriviamo 

(2.2) F(2)=) f_2" 

la condizione 2.1 diventa 

(2.3) F(Az) - AF(z) =0 > f(x" - 2) =0 


Si hanno dunque due possibilità: 
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n 
a) X # XA per wnl allora È =0 n> le F(z) assume la forma normale non 


risonante 
(2.4) F(z) = fe 


b) aa per n = gq + 1 con LeZ allora fo per n # £q + 1 e F(z) assu- 


me la forma normale risonante 
(2.5) F(z) = 25, + w(23)1 


dove w{x) è olomorfa in un disco e w(0) = 0. 
Notiamo che se |xl # 1 siamo sempre nel caso a) mentre se |\| = 1 
iw _. e” è 
allora detto XA = e siamo nel caso b) se >» è irrazionale, mentre se 
LP 


a con p,q primi siamo nel caso b). In questo ultimo caso il gruppo delle 


trasformazioni R(z) = Az è discreto e formato solo daq elementi. 


Coniugazione Diamo una mappa olomorfa non lineare vale il seguen- 


te teorema sulla coniugazione formale. 
Teorema. Sia F(z) una mappa olomorfa 
(2.6) z' = F(z) = Az + f(z) 


n : ; : ? 
allora se X = A per n > 1 risulta coniugata da una serie formale con la sua 


si 


" at è ; 
parte lineare, se \' “ = A allora è coniugata con la forma normale risonante 


espressa da 2.5. 
Dimostrazione. Sia 


(2.7) z = s(E) = E + é(E) 6(0) = 0'(0) =0 
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la trasformazione e &' = W(E) 1a forma normale. 


L'equazione di coniugazione 


(2.8) Fo®=® 07 
diventa nel caso in cui W(&) = E 
(2.9) (AE) - A@(E) = f(E + g(£)) 


Se allora scriviamo $(&) = B dn e" si ha 


nz2 
n - 
(2.10) (A - 2) dn a 
dove 
n 
(2.11) t,* a L ) dh + % 
s=2 kt ++k_=n 1 S 
k.21 
j 


e perciò i Tn dipendono solo da termini noti all'ordine n. Allora se \° # \ 


per n>1 l'equazione 2.1 determina dn Per nz2. 


Se invece "i = X per n = 2é + 1 chiaramente 2.10 non ammette solu- 


zione. Se però 
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(2.12) (E) = n& + E w(el) 


dove w(z) = E wz" allora l'equazione di coniugazione si scrive 
nz1 
o(1E)-A6(E) + em(e®) = f(e46(6)) - 
(2.13) 
- [olne + Ew(E )) - 6(AE)] 


Infatti da questa equazione si ha 


(2.14) (A 


o 


- + S 
1)è “ dn, 49+1 n 


dove i d dipendono solo da quantità note all'ordine n. L'equazione 2.14 è ri- 


solta da 


(2.15) W 


> -l 
ona * n° “ “i Sei 


Per quanto riguarda la convergenza valgono i seguenti risultati [5], [6] 


Teorema (Poincaré Dulec). Se |x|®1 allora $ è un diffeomorfismo 


analitico. 


Dimostrazione. Si basa sul metodo delle serie maggioranti. 0s- 
serviamo che essendo f(z) è olomorfa in un disco possiamo sempre operare una 
trasformazione in modo che sia UAN Detta quindi c una costante tale che 
flag” < c per n>2 e prendendo î moduli nelle equazioni 2.10, 2.11 si trova 
che (CA è maggiorata da c volte una somma espressa dal lato destro di 2.11 ove 


t, è sostituita da 1 e i da 4 . Ne segue che NS On dove 
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n=2 
2 
(2.17) ole) -E=c) o") =c Canto. 


che è risolto da 


— 
+e- Vite) -4€ (1+c) 
[aan otte) «SS 


e si vede che o(&) è olomorfa in un disco di raggio R = 1+2c - V1420)î-1. 


Teorema (Siegel [7]). Se |x|= 1 e x" = 1 per nd1 cioè rase!” 


con w/2r irrazionale e se 
n_,,-1 u 
(2.19) |a] <y n uzl 


allora esiste un diffeomorfismo analitico ®(&) che coniuga F(z) con la sua 
parte lineare. 
La dimostrazione di Siegel è ottenuta con il metodo delle serie 


maggioranti, noi daremo invece la dimostrazione usando il metodo di Newton. 


Teorema (Camancho [8]). Se |\] = 1e aa per n = 2q+1 cioè 
x = e'” con re = : » PsQ primi tra loro allora esiste un omeomorfismo locale 


® che preserva l'origine (9(0) = 0) tale che 


1 


®° o0Fo (E) = Az(1+ 259) 


per un K>1 intero. 
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Non c'è dunque coniugazione analitica ma solo un risultato di cor- 
rispondenza continua con le orbite di una particolare forma normale. 


Topologia delle orbite. Quando X = el è non risonante le orbite 


della forma normale sono curve chiuse 


È È 


Se invece w = 2 î le orbite della forma lineare sono date da q punti su un 
cerchio mentre quelle della forma normale sono date da curve emananti dall'o- 


rigine con 2q foglie 


Tale risultato dovuto a Julia può essere dimostrato costruendo i] 


flusso interpolante. Infatti se &' = W(E) = AE (1+ 59) si trova che in un si- 
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x 


stema rotante Z = e il flusso Z(t) ha un campo di vettori autonomo Z = U(Z) 


dove U(Z) = 284,1 + 0(Z)] e l'equazione dell'orbita che in coordinate polari 


è 


d 
» = p cotg gg 


ha per soluzione p = |sin qo|}/9 " 


3. IL METODO DI NEWTON 
Consideriamo l'equazione di coniugazione 
(3.1) W(AE) - Au(E) = f(E+U(E)) 
con |\|] = 1 e tale che 2.19 sia soddisfatta. 
Per risolvere la 2.19 usiamo un metodo di Newton modificato [4], [9]. 
Cominciamo col restituire 3.1 con l'equazione linearizzata 
(3.2) (AE) - Ag(E) = f(£) 
e definiamo ®(E&) = E+4(E) con g(0) = #'(0) = 0. Chiaramente 1a è non coniuga 
più la F con ja sua parte lineare ma con F_ (e) =\E + f1(8) dove f (e) è un re- 
sto. Essendo 
(3.3) ® Fog F 


si trova che f,(8) soddisfa l'equazione funzionale 


(3.4) f, (e) = AG(E) -o(xE + f1(8)) + f(E +9 (E)) 
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che usando la 3.2 diventa 


(3.5) f(e) = (AE) -@o(AE + f1(8)) + f(e+t0(2))-f(£) 


Il procedimento viene quindi iterato secondo îo schema 


In questo modo si ha una successione LITEZZZLA di trasformazioni 


e di nuove mappe Pi tali che detta A la trasformazione composta zh}? si ha 


ni L = iS 
(3.6) Yo F' o A Fn tn 0% ®, 0...0 d 


Ciascuna funzione LA è analitica in un disco P. dove tro} una successione mo 


notona decrescente con limite r_ > 0. Si prova che A SUP, è olomorfa e che 
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il limite per n>e è olomorfo in |z| < r_i inoltre la successione dei resti 
ui; converge-a zero per cui la successione Fi converge a Xz. 


Osserviamo che 


n 
(3.7) f_(2) = 012° *1) 


come ci si deve aspettare per i resti di un metodo di Newton. Infatti la rela- 


N4 
zione è vera per n=1 e supposto vera per n si trova $n4102) = 0(2° L; ed il 
risultato segue da 3.5 ove faf]3d son sostituiti da tt fn+1? Ènt1 poiché 
E i] 
(3.8) f+1(2) = fal2) è (2) [1+ 0(2)] 
n+tl 


Notiamo che la differenza con il metodo di Newton usuale è dello stesso ordi- 


ne del resto. Infatti riscrivendo 3.1 nella forma #(y) = 0 dove 
(3.9) F(Y1)=YoLl- Foy 


essendo L(&) =ME l'applicazione lineare. Allora detta I l'applicazione iden- 


tità e posto Y = I + y sviluppiamo il funzionale rispetto ad I univoco 
F(I + 4) =F(I) + F'(I)y += L-F' + gol - F'(1)y +2 = 


= -f + vol - Ly - [f'y +9] 
Dunque l'equazione linearizzata è ottenuta annullando vol - Ly - f. Osservando 
però che allora f'yw risulta essere un o(IvI?) come il resto 4 è chiaro che non 
c'è differenza rispetto al metodo di Newton stesso. 


Consideriamo allora l'equazione 3.1 che riscriviamo 


XXI-18. 
Osserviamo che M è un operatore lineare definito nello spazio A, 
delle funzioni olomorfe in |z{ < R, che è a un so1 valore sotto l'ipotesi 
ar per X>1. Se inoltre la successione 


(311) 


è superiormente limitata allora applica A, india con p* < p. Noi ci mettere- 


mo nell'ipotesi di Siegel 2.19 con u = 1 cioè 
(3.12) 1a - 1] <m 


questo al solo fine di rendere immediate le stime. Scegliamo in ad le norme 


u 
rs 
> 
3 
3 
AN 
eci 


(3.13) Ifl, 


e definiamo 


(3.14) s(r) = deh. DL ali. e 
r n 
nz2 
Allora si ha che è = Miro ta poiché 
Hi n n-1 
(3.15) lol, sii dear Peayvsenkr fil ‘yrés 
|a al 


Inoltre in un disco di raggio r(1-0) stimiamo la derivata 
n-l n-1_ y6 
= - - £ 
(3.16) I#l,(1-6) yInr If (n-1)(1-6) S 
PSE n -l 
poiché n(1-8) < (e6) per nz1. 


Per stimare il resto occorre premettere un lemma tecnico la cui 


prova è data in (10) 


XXI-19. 
Lemma. Se diveg e M = max( él,» Ivi.) ef <L.M allora 
1 ea ' a 
(3.17) IF(1+9) -F(140)], 5 1f',,y 10-41, 
Possiamo cioè enunciare il teorema sulla stima del resto. 


Teorema. Se fed, allora ti) dove 


(3.18) o = vY6 , mat 

l-e 
essendo e il numero di Nepero . Nel disco di raggio "i dato da 
(3.19) rg * r(1-0-n-6) 


dove o è arbitrario perché o + n +. o < 1 si ha la stima 


(3.20) If], se=-> 
dai 1 26 


Dimostrazione. Definiamo una sfera 4 nello spazio di Banach delle 


funzioni h SI tto) data da IR (1-5-n) s Tn , e consideriamo il funziona- 
le 


(3.20) Th(&) = o(ac) - d(ac + h(£)) + F(E+9(E)) - F(e) 
Mostriamo che T applica 2 in u ed è ivi contraente. Infatti usando 3.17 


[Th] s [o'| h + 
r(1-0-y) P(1-0-n) 4h (1-5)! Latiosal 


+ PI rn) (1.525) le (1-5-n) 
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NA 


Poiché Ihl L(1-0-n) rn e I (1-0-n) < lol, < y6r = ro 


usando 3.16 con 9 + o e la definizione di n data da 3.18 si ha 


ITAL, (1-0-n) fe 


D'altra parte se he® . he 


1 2le(1-0-n) s lo loti-a-n+ rn Ih,-h2l(1-0-h) 3 


IS In sE i 
Sto Irtaletizo-m Fe lrtalr(1-0-n) 


Schema iterativo. Usando il precedente risultato troviamo che lo 


schema iterativo è 


(3.21) If I, S $n+1 » Ea” re a) 
nt1 
dove 
x n 5 °n 
(3.22) Sn+t1 ecc gig Spa = i, 
2 n l-e 


con la condizione su 6 che ta; +0<1 


DI 
PE 2 Y n n RR , n ste 
Poiché di —— se si sceglie n ti si trova 8n+1 fi 


che è uno schema detto superconvergente. Noi sceglieremo invece 


(3.23) 5) =x 6 
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(3.24) avra ag La EI 
dove 
yYS ur 0) 
(3.25) Ga 1 * 
(o) 0) o Ca 1 (o) 26; ue Ly 


= <R. 
avendo posto I sr) dove la 


Si vede che . + o per n + © e si può dare una stima per r, con me 


lim pig A [I (1-9,7%n7%) > 
nro n=0 


Co] 


zr, Pf È logr1-x"{o, Paga DHr* 


n=0 
o © nX 

= ro expi -p pa È (ot * 6)% 2 

n=0 k=1 

PN Li 
CO (o +nto ) cuni 

al LR ail ti 

= r_ exp{- a L k } roi! d *W 90) 


Per ottimizzare la stima dobbiamo cercare il massimo di 
1 


A = lex _. 
A rl! 900 90) rispetto ad rv dove O<r fR e 0<x<1. 


Convergenza di n . La convergenza di questa successione si prova 


osservando che 
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N 
HW 
" 
> 
mm 
_ 
N 
mm 
_ 
N 
m 
+ 
> 
mm 
__ 
N 
aS) 
_ 


sicché si ha 


N 
" 


+ = + 

Zi + YnlZn) Fn tal) + + den) 
Ora notiamo che LA è olomorfa in |z| < fn-1 € la sua immagine è in 

un disco di raggio Yn-1 


lo 1 =I+6.1, sr + ol Srtr 0. Sr 
n ni, n n Fn=1(1-0 n-1 n-1 n-1 
Allora le composizioni 9-1 (o) LA sono ben definite e infine Unl2) converge per 


n + poiché 


n n 
bl, s) dal, s Vl 
me e o a alta) 


n n 
ki Me 
I rog "o di Lx si” SA lag 


Possiamo osservare che le stime ottenute sono anche quantitativa- 


Vi 
È 


mente discrete; infatti per il numero d'oro w = 2 si trova r, > 


81 
0.048 che va confrontato con il valore o = ..326 del raggio di Siegel otte- 
nuto con altri metodi [10], [11]. 


Estensioni. Il metodo qui descritto è stato esteso da S. Mermi al 
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[12], [14] al problema di Siegre in ci. 
Introducendo quindi opportune norme pesate e modificando la con- 
dizione diofantina è possibile una estensione anche al caso infinito di- 


mensionale d + ©, 


Ringraziamenti. Sono estremamente grato al Dr. Mermi per i suggerimenti e le 
osservazioni fatte su questo manoscritto e su tutte le informazioni bibliogra- 
fiche riguardo al capitolo 2 contenute in un suo seminario tenuto all'Universi- 
tà di Dublino [15]. 
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APPENDICE 


Diamo qui un cenno alla costruzione del flusso interpolante nel pro 
blema di Siegel risonante. Sia data la mappa 


Si, | ast (A-1) 


z' = T(z2) = e'[z+z 
con p,q interi primi tra loro. Consideriamo quindi un flusso Z = y(z.,t) con 
y(z,0) = z, interpolante l'orbita di T cioè tale che y(z,n) = 1" (2) yn in- 


tero. Sia inoltre V(Z) il campo vettoriale associato al flusso 


a |a 
ct IN 


= V(Z) (A-2) 


Si determina allora una soluzione a w(z,t) e V(Z) per via perturbativa; si scri- 
ve cioè 


sizat = rata. (A-3) 


vizy = nes bl Pa... (A-4) 
imponendo che (A-2) sia soddisfatta con le condizioni 
v(0,z) = z y(1,z) = T(2) (A-4) 


La struttura di gruppo del flusso y(t,z) assicura che l'interpolazione a t=1 


implica l'interpolazione per t = n intero. Sostituendo A3 in A4 si trova 


i ji(q+ +1 
V= iw el, + [iwa(t) + b(t)e!(0 l)ut Sa aiall PPP 


] 


, BL 
e imponendo che A-2 sia soddisfatta all'ordine 29° si ha 
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ji (q+ 

Î alt) = iw alt) + b(t) e'(9*Mot (A-5) 

da cui 
iwt t iqws 
alt) = e il b(s) e IS ds (A-6) 
Co) 
i % iqws 

Poiché (A-4) implica a(1) = e'’ si deve avere che f b(s) e dSds=1 e una 
scelta possibile è Ù 

b(t) = e '9t (A-7) 
Ne segue che 

y(2,t) = gut [z+t Pa Bo (A-8) 
e quindi 

MEO) do a I, (A-9) 


+ 
Un calcolo diretto di T"(2) all'ordine 2° * mostra che (A-8) interpola corret- 
tamente l'orbita della mappa. 
Se ci mettiamo nel sistema rotante definito da 
iwt 


Z=e € (A-10) 


il campo vettoriale W(E) associato al flusso discreto autonomo cioè 


dove 


XXI- 26. 


u(e) = e ®t(iuzev) = 8811 


Di sar 
A tutti gli ordini può essere fatta una scelta analoga alla (B-7) in modo che 
nel sistema rotante il flusso corrisponda ad un campo autonomo. 


Le orbite sono allora facilmente calcolabili poiché in coordinate 


polari le componenti di W sono 


+ + 
LA = o î cosqe , LA = p A sin gg e l'equazione dell'orbita è 


d 
FA = p cotg qe 


XXI-27, 
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